Low-frequency noise in tunneling through a single spin 
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We propose measurements of low-frequency noise in the tunneling current through a single 
molecule with a spin as an experimental probe for identifying a mechanism of the spin-dependent 
tunneling. A specific tail near the zero frequency in the noise spectrum is predicted; the amplitude 
and the width of being of the same order of magnitude as the recently reported peak in the noise 
spectrum at the spin Larmor frequency. The ratio of the spectrum amplitudes at zero- and Larmor 
frequencies is shown to be a convenient tool for testing theoretical predictions. 
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Tunneling currents via a microscopic system, such as a 
quantum dot, or a molecule or an atom with a localized 
spin 0,1^1^ attracts considerable attention in the context 
of the problem of quantum information processing. The 
tunneling current depends on spin-dynamics and thus en- 
codes its features; at the same time, the tunneling (mea- 
suring) current influences the spin dynamics itself. Thus 
the tunnehng via a single spin current measurement can 
provide information on spin orientation and its dynamics 
and offer an example of an indirect-continuous quantum 
measurement 0. 

Scanning tunneling microscopy (STM) experiments 
on a single molecule with a spin in the presence of a mag- 
netic field Bq, have revealed a peak in the current noise 
power spectrum (i. e., in the current autocorrelation func- 
tion) Viuj) at the Larmor frequency lul = jBq, where 7 is 
the gyromagnetic ratio. Experiments were performed at 
room temperatures and found the signal-to-noise ratio R 
(the ratio of the power at peak frequency to the shot noise 
power) to exceed unity and be almost independent of the 
orientation of the applied magnetic field Bq. Experimen- 
tal results are hard to explain within the framework 
of a single-spin non-rclativistic model since electrons in 
the leads are polarized by the same magnetic field that 
acts on the spin i.e. along Bq, and do not couple with 
the oscillatory, components of the spin, which arc per- 
pendicular to Bq. a possible relevance of the spin-orbit 
interaction has been suggested 0] ■ Recently, Levitov and 
Rashba noticed that in the systems with the low space 
symmetry (such as a dot or a molecule near the surface) 
the nonvanishing orbital moment of tunneling electrons 
couples them to the mediating spin. They suggested that 
this mechanism may lead to a significant effect of the spin 
oscillatory component on the tunneling current. 

Yet to understand experimental results not only the 
existence of the peak in the current power spectrum is to 
be explained, but also the origin of the large, R > 1, 
signal-to-noise ratio and its weak dependence on the ori- 
entation of Bq . An important step towards formulation 
such a model has been done by Bulaevskii, Hruska, and 



Ortiz who included in the model the non-rclativistic ex- 
change coupling of a single spin 1/2 and the tunneling 
electrons 0, 0|- Their approach that followed was 
based on the Keldysh formalism and the Majorana- 
fermion representation |^ for the spin, thus taking into 
account the nonequilibrium effects in spin dynamics ex- 
plicitly. They found the spin distribution function and 
the current-current correlation function and discussed 
the dependence of R and line width F in the current 
power spectrum on the applied voltage V between leads, 
the applied magnetic field Bq, the temperature T. They 
also obtained R- and F-dependence on the degree and ori- 
entation of electron polarization in the right- {a = 
R) and left (a = L) current-leads in the steady state (this 
state establishes during the transient time after the volt- 
age or tunneling matrix elements are switched on) 0] . 
The results of |3 explained several qualitative features 
of both average tunneling current through the spin and 
noise spectrum at the Larmor frequency. The quantita- 
tive agreement with the experiment was not achieved, yet 
the model of is attractive and warrants to be explored 
further. The task now is to identify experimentally ac- 
cessible effects that could test the underlying physics of 
the tunneling-through-a-spin phenomenon. 

We propose measurements of the low-frequency noise 
(LFN) in the tunneling current as such a probe. In this 
Letter we develop a theory of LFN of the tunneling cur- 
rent adapting the Ref. model. We predict a tail near 
zero frequency in the noise spectrum having the width 
of the same order of magnitude as that of the peak at 
lul. The LFN is expressed through the same quantities 
as the noise at the the Larmor frequency, and thus the 
ratio p = V{0)/'P{ujl) turns out to be a function of the 
bias voltage V, magnetic field Bq, as well as of polar- 
ization of the leads and of tunneling coupling. Thus the 
experimental study of behavior of the parameter p offers 
a unique tool to check on our understanding of tunneling 
through a localized spin. 

Low-frequency noise of the tunneling current. We use 
the same notations as in 0: voltage is measured in the 
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energy units, thus we write just V instead of eV, further- 
more B stands for gfisB, T stands for ksT, and lo repre- 
sents hio. Thus in our notations B = cj^. The Hamilto- 
nian is that of the two- leads Kondo model P, |3i 13 where 
the direct tunneling term is also included: 



Ti. — He -\-'Hs +'Hti "Ht — Ti-ref + Ti-tr, 
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where cj,„g. (Cq,„^) creates (annihilates) an electron in the 
left or right lead (depending on a e {L, R}) in the cigcn- 
state n, and with spin a. Further, = e„ — where 
e„ is the energy in the state n and is the chemical po- 
tential in the lead a, while 3 represents the three Pauli 
matrices. ^^L^ ^rr and tI^^^ are tunneling matrix el- 
ements due to the exchange interaction for the electron 
tunneling from the leads to the molecule with the spin 
1/2, while To is the direct tunneling matrix clement. We 
take them as real numbers. The spin localized in the 
molecule is described by the operator S = {Sx, Sy, Sz)- 
Figurensketches the physical setup we want to study and 
which basically represents the model Hamiltonian Ti.. 
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FIG. 1: Schematics of the physical system. The electronic 
tunneling current is established by a dc voltage V. 

The electrical current operator can be written as 
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clnAt)iTRL)aa'C^n'.'{t)+H.C.. (3) 



Since the spin-dependent tunneling amplitude Tul, 
Eq. 10); contains two terms, the current can be schemat- 
ically presented by two vertices: Tq = Tq Sa-a-' corre- 
sponding to the spin-independent tunneling, and Ts = 

(ex) 

T^L S • (Ta-a' corresponding to the spin-dependent part. 

(ex) 

In the following we assume Trl ^ ^o- Since we are 
interested only in low frequencies, lu <^ B,V , only z- 
component should be kept. Thus the current noise can 
be expressed by the schematic diagram shown in Fig. |21 
The left and right electron blocks arc nothing but the 




FIG. 2: Skeleton diagrams for calculation of the low-frequency 
noise. The wavy line corresponds to the correlation function 
{5Sl)^. 



derivatives dI/d{Sz), while the wavy line corresponds to 
the correlation function {5Sl)^ = {[5Sz{t),5Sz{Q)\ + )u]- 
This way we arrive at the natural expression for the low- 
frequency noise, 

V{uj) = {[Sm, <5/(0)] + ). = {dI/d{Sz)f {5Sl)^ . (4) 

To be concrete let us restrict ourselves to the case of 
fully polarized electrons in the leads. Then there exists 
the average current proportional to the average spin, (S), 
and given by the expression 

I{V) = /o(^) + Is(^)-(S), (5) 
/o(y) = ^e(l + niR • niL) T^pl V, (6) 

ls{V) = 2^e(mR + niL) To^lVo ■l^- (7) 

In the above equations po is the density of states per 
spin (DOS) of the leads at the Fermi level (when leads 
are different = PgP^ where is the DOS in the lead 
a); xo-a means the direction of electron spin polarization 
in the coordinate system with z-axis parallel to the total 
magnetic field, B = Bq + By. Here Bq is the external 
magnetic field while B^ is the additional dc magnetic 
field produced by tunneling electrons. 

The derivative dI/d{Sz) = Isz{V)^ see Eq. (O, may 
be found in Q, and we turn to calculation of the spin 
correlation function. In the equilibrium, the fluctuation- 
dissipation theorem yields: 



{5St)^^coth—lTnx{u:) 
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Here {Sz)t is the equilibrium average spin, F, is the de- 
cay rate for the Sz fluctuations. 

Far from the equilibrium the fluctuation-dissipation 
theorem generally speaking does not apply. However, 
in a non-equilibrium, but yet stationary state one can 
still use the result of formula © with the appropriate 
expression for the average spin, (Sz) = (l/2)/i/(6). 
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Here hf{h,v) is the function calculated in 0]: 

26(1 - mR,TOLz)-2u(TOR, - TOLz)+fe 01 



hf = - 



0+(l - mR^TOLz)-^ (rriB^z - mi^z) + 4>ib) 
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(t){b) = 6coth(6/2), 0=* 



{v + b) ± (j){v - b) . (10) 
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Equation H10|l makes sense provided m^j. ^ 1 and 
TTiLz 7^ il- Otherwise, if torz = 1, ttilz = ±1, the self- 
consistency equation is an identity and the spin steady 
state can be any. 

Equation can be derived, e. g., using the the tech- 
nique developed by Abrikosov , where a 1 /2-spin was 
interpreted as a pseudo-Fermion with the Green function 
g±\e) = {e ^ B /2 — X + iS)~^ . Here A is an auxiliary 
"chemical potential" which is send to infinity eventually. 
This trick allows one to remove extra unphysical states 
that appear because Fermi operators have more extended 
phase space than spin operators. The method was elab- 
orated by Maleev for the case of dynamical defects in 
glasses Schematic diagrammatic representation of 

the correlation function {6S'^)lj is given in Fig.Bl Deriva- 
tions can be carried out similarly to those in [13 where 
electron dephasing rates due to pseudo-spin defects were 
calculated. A similar procedure was also used for calcu- 
lation of the energy relaxation time of the electrons in a 
thin wire due to magnetic impurities . The correlation 
function is given by the diagram shown in Fig. |31a, where 
thick dashed line corresponds to the auxiliary fermion 
(7±(e). An important feature of the calculation is that 
for calculation of {SSD^ one has to calculate the self 

(ex) 

energy up to the fourth order in Tpr (Tt- as shown in 
Figs. and , see also Appendix in . The wavy 
line, Fig.|3|i, corresponds to Sx — Sx propagator and solid 
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FIG. 3: Schematic diagrams for calculation of {5S'^)uj- Dashed 
lines represent auxiliary fermions, while the vertices Sz and 
Sx represent the tunneling amplitudes T'^^az and T^^l^ax, 
respectively. 



lines represent tunneliiig electrons. This propagator was 
actually calculated in [3. A similar complication arises 
when calculating the Sz (namely, T^^azSz) vertex part, 
see Fig. |3J;,f. When niR = niL (except for maz = ±1), 
the result is the same as for unpolarized electrons 0: 



hf{b, v) = tanh | — 



^{2 + 0) 



(11) 



Here 6 = (T^^ + Tl^)/Tl^- We would like to empha- 
size that matrix elements T^a describe electron tunneling 
from the lead a to the spin site and back while T^r de- 
scribes spin mediated tunneling between the leads. Thus 
the value of 9 is extremely sensitive to the location of the 
spin with respect to the leads. For the perfectly symmet- 
ric configuration Tll ~ Tj^n — Tlr, therefore 6 = 1/2. 
Since both Tll and Thr increase dramatically with a de- 
crease of the distance from the spin to the corresponding 
lead so does 9 with an increase in asymmetry. As a result, 
in the asymmetric configurations the average spin is ac- 
tually equal to its equilibrium value. This was noticed, in 
particular, for the similar problem of the electron tunnel- 
ing mediated by the presence of the structural two-level 
system in . However for symmetric configuration the 
spin coupling to the electrons tunneling between the leads 
is as strong as its coupling to the electrons in any of the 
leads. As a result, the average spin is controlled by the 
combination of electron energy distributions within both 
of the leads, and it is out of equilibrium provided u > 1. 
The tunneling electrons reduce the spin magnetization 
which drops as 1/V at large V. 

Substituting Eq. into Eq. ^ we obtain 



2r, 



/ dl 



^2 + r2 \d{Sz 

Using the expression from 
of the current. 



[l~h}iv,b)] 



(12) 



for the spin-dependent part 



dI/d{S) EE h{V) = 2^e(mfl + mL)roT^°^"Vo^ , 

we arrive at the final expression for the low-frequency 
noise: 

2r, 



Pov' 1 



hjiv,b)] 
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Po 



{muz + rriLzY 



(13) 
(14) 



This result agrees with calculations by Shnirman et 
al. 01 . According to Eq. ifT^ . the effect is strongly de- 
pendent on a degree of electron spin polarization. Here 
we note that the polarization can arise not only due 
to spin-dependent tunneling amplitude, but also as a 
result of electron motion through the molecule where 
the localized spin is located. Indeed, as it is known 
for semiconductor structures 0|, the electrons tunnel- 
ing through the barriers with no inversion center be- 
come spin-polarized if the tunneling electron has a com- 
ponent of wave vector fcn parallel to the barrier plane. 
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The difference between the tunnchng exponents for oppo- 
site spin direction can be estimated as ~ ^{2mk\\lh^)kd 
where 7 characterizes an efficiency of spin-orbital inter- 
action (which for typical semiconductors is of the order 
of 10"'^^ erg- cm'^) while d is the tunneling length. Of 
course, there is a difference between the organic molecule 
and a semiconductor. However the apparent absence of 
any pronounced symmetry in the considered situation 
can result in a finite k\\ , which will, in turn, lead to a 
spin polarization. In particular, for fc|| ~ 10^ cm""'^ (two 
orders of magnitude less than typical atomic value) and 
for the values of 7 two order of magnitude less than the 
estimate given above the degree of spin polarization can 
be of the order of 10~^ which is larger than the spin polar- 
ization of electrons in metals at i? ~ 1 T. Thus the factor 
mentioned above can significantly increase the coupling 
of the tunneling electrons with localized spin. 

Comparison with the noise at Larmor frequency. Ac- 
cording to Eq. H13(l . the low- frequency noise is repre- 
sented by a Lorcntzian tail with the width F^. At w — > 



as parameter b increases. Thus, the ratio p can be of the 
order of unity. 

In conclusion, we have calculated the low-frequency 
noise power in the tunneling current, 'P(O), and demon- 
strated that the ratio p{V,T) = P(0)/7'(wl) is the uni- 
versal function of dimensionless voltage and magnetic 
field, eV/kBT and g^BB/ksT, respectively. This opens 
the route for identification underlying mechanisms of the 
noise in tunneling current by comparison of the mea- 
sured dependence of p(y, T) upon voltage and tempera- 
ture with the obtained J-{V, T). We also noted that spin 
polarization can be obtained due to intrinsic spin polar- 
ization while tunneling through a complex molecule. 
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pointed out some mistakes in our original draft and made 
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Let us compare this result with the maximal value of 
the noise at Larmor frequency, Vl- According to 0, 
the magnitude of the noise near the Larmor frequency at 
> -B is given by the expression 



Pu,=Pl[v^+vbhf{v,b)] 



Ti + iuj- By 



(16) 



Pi = (^rorsp^r)^imfli + mLxr 

(for fully polarized electrons). At V < B , Pu, 
Comparing now noise magnitudes at a; - 
Lu — > lul respectively at V > B, we have 



= 0. 

and at 



— / LiJ\UJ—LJl 



Po_ 
Vl 



{Pi/r^)[v^+vbhf{v,b)], (17) 

2 

T{v,b), (18) 



rriLz 



T{v, b) = [l- h){v, b)]/[l + {b/v) hfiv, b)] . (19) 

It rather difficult to provide realistic estimates for the 
ratio Tj_/Tz since spin relaxation and dephasing can be 
produced both by tunneling electrons and by some de- 
grees of freedom in the leads. The contributions of the 
tunneling electrons to r_L and are calculated in the 
Bloch-Redfield approximation in Ref. What is impor- 
tant that in general Fj^ > Tz/2, and the ratio F^/F^ can 
be measured experimentally. The dimensionless function 
J-{v, b) grows monotonously with v for a given b, it tends 
to 1 at w &, the plot J-{v) shifts slightly downwards 
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